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In recent years, Loop Quantum Gravity has emerged as a solid candidate for a non- 
perturbative quantum theory of General Relativity. It is a background independent 

o. 

theory based on a description of the gravitational field in terms of holonomies and 

fluxes. In order to discuss its physical implications, a lot of attention has been paid 

OS ' to the application of the quantization techniques of Loop Quantum Gravity to sym- 

CN 

metry reduced models with cosmological solutions, a line of research that has been 

O . 

O -1 called Loop Quantum Cosmology. We summarize its fundamentals and the main 

(3jT)! differences with respect to the more conventional quantization approaches employed 
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in cosmology until now. In addition, we comment on the most important results that 
have been obtained in Loop Quantum Cosmology by analyzing simple homogeneous 
and isotropic models. These results include the resolution of the classical big-bang 
singularity, which is replaced by a quantum bounce. 

PACS numbers: 04.60.Pp,04.60.Kz,98.80.Qc 



1. MOTIVATION 

Gravity is the only fundamental physical interaction which is not yet satisfactorily de- 
scribed quantum mechanically. Even without adhering to the belief that all fundamental 
interactions should finally be unified in a single theory, a strong motivation to search 
for a quantum theory of gravity comes from the very own results of General Relativity. 
The classical singularity theorems that arise in Einstein theory [l| imply that (in a vari- 
ety of physically relevant situations) the predictability breaks down, so that the regime 
of applicability of General Relativity has been surpassed. Therefore, a new and more 
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fundamental theory is needed for a correct physical description. 

In trying to quantize General Relativity, the first obstacle that one finds is that Einstein 
theory is not renormalizable as a quantum field theory, so that a conventional perturbative 
quantization cannot be performed. In this context, an alternative quantization program, 
joop Quantum Gravity (LQG), has recently been proposed for General Rel- 
LQG is an attempt to construct a nonperturbative quantum theory of 
gravity using techniques similar to those of gauge field theories (e.g., Yang-Mills). The ap- 
plication of these nonperturbative quantization techniques to simple gravitational models 
with application in cosmology, such as homogeneous and isotropic spacetimes with differ- 
ent types of matter content, has given rise to a new branch of gravitational physics called 
Loop Quantum Cosmology (LQC) |5|]. 

2. HAMILTON! AN FORMULATION OF GENERAL RELATIVITY AND 

ASHTEKAR VARIABLES 

LQG is a nonperturbative canonical quantization of General Relativity; therefore, it is 
constructed starting from a Hamiltonian formulation of Einstein theory 6] . Let us review 
very briefly this formulation. 

We consider globally hyperbolic four- dimensional spacetimes (g a p, M = R X £), where 
g aj 3 is a Lorentzian metric, Greek indices are spacetime indices, and S is a three- 
dimensional manifold. For General Relativity, once E is given, the physically relevant 
information to determine the classical solutions is contained in the spatial three-metric 
h ab induced on S, and in the corresponding extrinsic curvature K ab = ^C n h a b, where 
Latin indices from the beginning of the alphabet denote spatial indices, n is the unit 
normal to E, and C is the Lie derivative. Equivalently, we can adopt co-triads e l a (rather 
than spatial metrics) to describe the system. This allows the coupling to fermionic matter 
fields. With 77^ being the Euclidean three-metric, we have the relations 

h ab = ei Vij 4, Kl = K ah e)^. (1) 

Here, the triad e b is the inverse of the co-triad, e l a e\ = 5 b a and e % a e°- = <5], and Latin 
indices from the middle of the alphabet are internal SU(2) indices, corresponding to the 
symmetries of the Euclidean metric rjij under linear transformations. 
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A set of canonical variables for General Relativity (in the sense that their Poisson 
bracket is proportional to the identity) is given then by the densitized triad Ef and the 
extrinsic curvature in triadic form Kl: 



a ' 



E:-.= Vtethel K - {K:(x),E b 3 (y)} = 5 b a 5^\x-y). (2) 

In this expression, x and y are two generic points in E, and we have chosen units such 
that 8ttG = 1, where G is Newton constant. 1 

Actually, we can replace the triadic extrinsic curvature by a connection valued 1-form 
taking values on su(2), and still obtain a canonical set of variables. For this, it suffices 
to realize that the co-triad determines an su(2)-connection compatible with it, T l a , and 
notice that the sum of this connection with any vector (both from the internal and spatial 
viewpoints) provides again an su(2)-connection valued 1-form. Therefore, at the classical 
level, we can simply replace K l a with A^' l a = T l a + 7-fQ. Here, 7 is a nonzero constant 
called the Immirzi parameter, and its presence can be seen to lead to an ambiguity in the 
quantization [?], Q] which is usually resolved in LQG by appealing to the recovery of the 
Bekenstein-Hawking law for the entropy of black holes 9|. For simplicity, we will set it 
equal to one from now on. The calculations for general 7 can be easily reproduced along 
the lines explained below. 

We will thus adopt as canonical variables the set formed by A l a = T l a + i^!and Ef. In 
General Relativity, these variables are subject to three types of constraints |0, |f|. First, 
there is a Gauss constraint which generates SU(2)-transformations, 

Qi := d a E? + e {j k AiE a k = 0. (3) 

In addition, the invariance of the theory under spatial diffeomorphisms is reflected in the 
so-called vector or diffeomorphism constraint, 

V a := F*fi = 0, (4) 

where F* b is the curvature of the connection A l a , namely 

F i ab = 2d [a A\ ] + e i jk Ai a At (5) 



1 In the following, we also set h~c^ 1. 
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Here, e^k is the totally antisymmetric symbol. Finally, the invariance of General Rela- 
tivity under time reparametrizations leads to a scalar constraint, also called Hamiltonian 
constraint, which in vacuo takes the expression 

S := E«E) (e Va* - 4*f.4) = 0. (6) 

Given the four-dimensional covariance of Einstein theory, General Relativity is a com- 
pletely constrained system, i.e., the total Hamiltonian which generates the dynamics is 
just a(n integrated) linear combination of constraints. In particular, apart from boundary 
terms, the Hamiltonian vanishes on classical solutions. On the other hand, it is worth 
pointing out that General Relativity is formulated in terms of connections and densitized 
triads without introducing any metric background structure. This background indepen- 
dence plays a fundamental role in the theory and will be a basic guideline for the selection 
of a quantization procedure in the construction of LQG. 



3. HOLONOMY AND FLUX ALGEBRA 



Since SU(2)-transformations are symmetries of our gravitational systems, only the 
gauge invariant information about the connection is physically relevant. Taking into 
account that this information is captured by the Wilson loops Q, we can then 
replace the connection by SU(2)-holonomies. More specifically, from now on we will 
consider holonomies along piecewise analytic 2 edges e, where we understand that an edge 
is an embedding of the interval [0,1] in E 0, [3]. We call h e the corresponding holonomy, 



h e = Pexp J A l a Tidx a . (7) 

Here, the symbol V denotes path ordering, and {jj = — \<Jj\ j = 1,2,3} is a basis in 
the algebra su(2), with aj being the Pauli matrices. Let us notice that the line integral 
appearing in the holonomies implies a one-dimensional smearing of the connection, and 
that no use of background structures has been made in the definition of the holonomy. 

Since the most relevant field divergences in our theory are expected to come from the 
appearance of the three-dimensional delta function in the basic Poisson brackets between 



2 



The restriction of piecewise analyticity ensures that the intersection between edges, as well as the 



intersection of an edge with a (piecewise analytic) surface, occurs in a finite number of points 12 1 
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our variables, and we have already smeared the connection over one dimension, it seems 
natural to try to smear now over two dimensions. Once again, we want to carry out 
this smearing without employing any background structure. Remarkably, this requirement 
can be fulfilled because Ef is a vector density. Hence, for any piecewise analytic surface 
S and any su(2)-valued smooth function /' on it, we introduce the associated flux of the 
densitized triad, 

E(S,f)= [ E?fe abc dx b dx c . (8) 

The defined holonomies and fluxes form an algebra under Poisson brackets. In the 
following, we take this algebra as our algebra of elementary phase space variables. From 
this perspective, the quantization of the system amounts to constructing a representa- 
tion of this algebra. A keystone result in LQG is a uniqueness representation theorem 
known as the LOST theorem (after the initials of its authors [3]). The LOST theorem 
states that there exists only one cyclic representation of the holonomy-flux algebra with 
a diffeomorphism-invariant state (interpretable as a "vacuum"). Therefore, the choice 
of the algebra of elementary variables, motivated by background independence, together 
with the identification of diffeomorphism invariance as a fundamental symmetry suffice 
to pick out a unique quantization (up to unitary equivalence). 

In order to gain insight into the kind of quantization adopted in LQG, let us first 
call cylindrical those complex functions of the connection that depend on it only via the 



holonomies along a finite number of edges (forming a graph |12l|). We can identify the 
commutative unital *-algebra of these functions as the algebra of configuration variables. 
By completing it with respect to the sup-norm 3 (i.e. the supremum norm), we obtain 
a commutative C*-algebra with identity. Gel'fand theory ensures then that this algebra 
is (isomorphic to) that of continuous functions on a certain compact space, A, which 



is usually called the spectrum [14j. Smooth connections are dense in this space A of 
quantum generalized connections. Besides, the Hilbert space of any representation of the 
C*-configuration algebra is of the form L 2 (A, li) for some measure 11. The LOST theorem 
guarantees that there is a unique Hilbert space L 2 (A,lial) supporting a representation 
not just of the holonomies, but of the whole holonomy-flux algebra, and such that lial 



3 The use of this norm is motivated by the fact that, in a representation in which configuration variables 
acted by multiplication, the operator norm would coincide with the sup-norm. 
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(the so-called the Ashtekar-Lewandowski measure) is a diffeomorphism-invariant, regular 
Borel measure. This representation turns out not to be continuous and, as an important 
consequence, the connection itself cannot be defined as an operator valued distribution 

0. 



4. LOOP QUANTUM COSMOLOGY: FLAT FRW MODEL 

LQC confronts the quantum analysis of cosmological systems by applying similar quan- 
tization techniques to those described for LQG. Here, we will focus our discussion on a 
simple but physically relevant model, namely, the case of homogeneous and isotropic flat 
(FRW) cosmologies. As the matter content, we will consider a massless minimally coupled 
scalar field. 

The spatial manifold £ is topologically M 3 , endowed with the action of the Euclidean 
group. One can introduce a fiducial flat co-triad, °e l a , with the corresponding fiducial 
metric and triad, °h a b and °e°. Given the non-compactness of S, we also choose a 
reference cell adapted to the fiducial triad in order to integrate homogenous quantities, 
such as the symplectic structure or the Hamiltonian, without introducing infinities in our 
formalism. We use the symbol Vq t° denote the fiducial volume of this cell. Actually, 
physical results can be proven independent of these choices under a suitable definition of 



the elementary variables 15l.ll6j|. In more detail, one can fix the gauge and diffeomorphism 
freedom so that 

Aa = cV -V3 o e i r . ? and E a =pV ~ 2/3 Vdet h ° e y. (9) 

Here, c and p are constant on E (but not under evolution), and describe the only remaining 
degrees of freedom in our basic variables. In the following, we call T the subspace of the 
gravitational phase space (for full General Relativity) defined in this way. 
The gravitational symplectic structure induced on £ is just 

Q s = 3dcAdp, so that { c ,p} = -. (10) 

3 

The variables c and p are hence canonical, apart from the factor of 1/3. Holonomies along 
straight edges /i°e" in the fiducial directions suffice to separate symmetric connections, 
i.e., given two different connections, there always exists an edge of this kind for which the 
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corresponding holonomies differ 14|. We thus restrict our attention to those holonomies, 
ho e .(fi), which have the form 

/ i o ei (/,) = cos(^)l + 2sin(^)r J . (11) 

Similarly, densitized triads can now be smeared just across squares with edges parallel to 
the fiducial directions, 

E(SJ)=p°A(SJ)V - 2/3 . (12) 

The factor °A(S, f) measures only the fiducial area of S weighted with an orientation 
factor. In this sense, fluxes are totally determined by p, which therefore plays the role 
of a momentum. The configuration algebra, on the other hand, is generated by sums of 
products of matrix elements of holonomies. Thus, it is the linear space of continuous and 
bounded complex functions in M provided by finite sums of the form f(c) = Y2 n fn e* MnC . 
Its completion with respect to the sup-norm is known to be (isomorphic to) the Bohr 
C*-algebra of almost periodic functions [ijj]. 

5. BOHR COMPACTIFICATION AND POLYMER REPRESENTATION 

As we have commented, the configuration C*-algebra is the algebra of almost periodic 
functions. The (Gel'fand) spectrum of this algebra is the Bohr compactification of the 
real line, IRsohr [14]. This compactification can be understood as the set of group homo- 
morphisms from the group M (with the sum) to the multiplicative group T of complex 
numbers with unit norm. So, every x G Msohr is a map x : M —>■ T which satisfies 

x(0) = 1, x(pi +p 2 ) = x(p 1 )x(p 2 ) Vpi,p 2 eM. (13) 

Since T is a commutative group, the operation xx(p) := x(p)x(p) provides a commutative 
group structure in Msohr- This group is compact with respect to the Tychonoff product 
topology. We recall that the Tychonoff topology is the weakest topology for which the 
functions F p : J R-Rn hr — > T given by evaluation at p [i.e. F p (x) := x(p)} are all continuous 



for any p G M 14J ]. Besides, the real line is actually dense in IRsohr- This result follows 



from the fact that the algebra of functions f(c) considered at the end of the previous 



section separates points c G M 14] 



The compact group IRsohr is equipped with a normalized invariant measure under the 
group operation, namely, the Haar measure fin- The representation of the holonomy-flux 



algebra for LQC is given precisely by the Hilbert space L 2 (MBohr, Hh)- In addition, since 
p H is invariant under multiplication in the group, we get that, \/x G Msohr-, 

[1 - x{p)\ I F p {x)dp H {x) = 0, (14) 

from where it follows that 

/ F p (x)dp H (x) = 5° p . (15) 

Taking into account that, from our definitions, F Pl F P2 = F Pl+P2 and F* = it is 

straightforward to conclude that the set {F p ,p G M} is orthonormal (hence, the Hilbert 
space L 2 (lR,Bohr, Hh) is nonseparable). One can also see that this set is dense [uj ]. As 
a consequence, the Hilbert space L 2 (lR Bohri p H ) is isomorphic to the so-called "polymer" 
space of functions of p G M that are square integrable with respect to the discrete measure. 
The isomorphism is given by I : F p — > \p) \/p G M. 
Employing then the orthonormal basis 

{\ P ); P eM,(p\p) = 6§}, (16) 

and introducing the notation := exp (i/xc/2), the polymer "momentum" representation 
is determined by the following action of the holonomy and flux operators: 

p \p) = | \p) , iV M \p) = \p + ii). (17) 

In this representation, states take the general form 

|V> := 1>(p) b> ; ^|^(p)| 2 <oo. (18) 
pen pen 

Note that normalizable states ip(p) differ from zero only on a countable subset of the real 
line for the label p, because the sequence {^(p)} is square summable. On the other hand, 
it is worth noticing that the representation (of N^) fails to be continuous in /i, as can be 
easily seen by realizing that the state \p) is always orthogonal to \p + p) regardless of the 
value of p, ^ 0. Therefore, the connection operator c is not well defined, in total parallelism 
with the situation discussed for LQG. The failure of continuity makes evident that the 
representation is inequivalent to the standard Schrdinger one of geometrodynamics 
(often called the Wheeler-DeWitt representation). This lack of continuity explains why 
the Stone- von Neumann uniqueness theorem of Quantum Mechanics does not apply [17J , 
allowing the results of LQC to differ radically from those -physically unsatisfactory- 
attained in geometrodynamics. 
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6. QUANTUM FRW MODEL 



With our symmetry reduction to the flat FRW model and our choice of fiducial struc- 
tures, the triad adopts the expression e™ = sign(p)|p| _1 / 2 V^ 3 °ef. This triad diverges at 
the big-bang singularity, corresponding to p = 0. In the quantum theory, on the other 



hand, p has just a point spectrum |18| which coincides with the whole real line, since the 
basis states \p) have unit norm Vp e M. Since zero is included in this point spectrum, 
the related (inverse) operator is not well defined. However, it is actually possible 
to define a triad operator in terms of our elementary ones [^(J. Classically, we have the 
following identity V/2 G M: 

= i tr fcr'hoM {<(/,), yfiplf) • (19) 

Here, ho £i is again the holonomy along the edge °ej, and the symbol tr denotes the trace. 
Then, replacing Poisson brackets with — i times commutators, we obtain 



6 



sign(p) 



N^p^N^-N^N^. (20) 



It is not difficult to check that this operator is diagonal in the p-basis. Furthermore, it 
is bounded from above, so that the classical divergence at p = disappears quantum 
mechanically with this regularization of the triad [20]. In fact, this triad operator is such 
that it annihilates the state \p = 0). 

Since we have already fixed the gauge and diffeomorphism freedom, the only constraint 
remaining in the system is the Hamiltonian one. For flat FRW spacetimes with a massless 
scalar field <t> (and unit lapse), this constraint can be obtained from the evaluation of the 

nn 

following expression in the symmetry reduced model [15|, [16( 



H:=\! |detiT 1/2 (Pi - e ij k E^F* b ) = 0. 



(21) 



To define the operator corresponding to |deti5| (or to |deti£| 1 ^ 2 e tJ k E^Ej in the gravi- 
tational part of the constraint 20|) we proceed as we have explained above when discussing 



the triad operator. In addition, to introduce an operator representation for the curvature 



F^ b , we first recall the classical relation 



^ = - 2 J i 2i tr (^5^ Tt ^ V 'l' (22) 



which is valid for any real value of /i and where 

:= ho e Xft)ho ej (p,)h^l(p,)h^l (p,). 



(23) 



Nonetheless, after substituting classical holonomies by their quantum counterparts, the 
limit of zero regulator fi cannot be taken in the resulting curvature operator. This cir- 
cumstance is interpreted as a manifestation of the fact that, in LQG, the area spectrum 



211 ] . so that the square with edges [i°ei 



is discrete with a minimum nonzero eigenvalue 4], 
and j2°ej, employed to define hu,j](jl), cannot be shrunk to zero. The regulator is then 
fixed by demanding that the physical area of this square equals the minimum nonvan- 
ishing eigenvalue allowed in LQG, which we call A from now on. Hence, one gets the 
operator relation fi 2 \p\ = A. 

At this stage, it is convenient to relabel the p-basis by introducing the affine parameter 
associated with the vector field |/2d p [lfj]. This vector field can be regarded as that 
corresponding to the exponent ~//c in the holonomy N^. Taking into account that the 
physical volume of the fiducial cell is given by the operator V = |p| 3//2 , the above relabeling 
leads to a basis of volume eigenstates \u), where v = 4 sign(p)|p| 3 / 2 /-\/A. The operator 
Nfi is then defined to produce a constant unit shift in the new label [161 ]. 



N^u) := \u + l). 



(24) 



Using the standard Schrdinger representation for the matter field, so that the total 
Hilbert space is the tensor product of the polymeric one and of L 2 (M,d(p), and adopt- 
ing a suitable factor ordering, one finally arrives at the following quantum Hamiltonian 
constraint: 

3/2 — - — 3/2 

(25) 



H 



1 



n 



1 



4VAi 



1 



-6ft 2 + Pi 

-V2 _ 



1 



IP I 



N 2fl - iV_ 2A ) sign(p) + sign(p) ( N 2il 



2f± 



-1/2 



(26) 



The symmetric factor ordering adopted for Q arises naturally from the consideration of 
homogeneous but anisotropic models of Bianchi I type, where the ordering is well moti- 
vated, regarding the flat FRW cosmolog 1GS clS db special case with vanishing anisotropies 
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221 ] . It is straightforward to check that the above quantum constraint annihilates the 
state \p = 0) (or equivalently \v = 0)) and leaves invariant its orthogonal complement. In 
the search for nontrivial solutions to the constraint, one can then restrict all considera- 
tions to this orthogonal complement, so that the classical singularity, corresponding to 
p = 0, can be removed from the kinematical (gravitational) Hilbert space 22j, |23J]. In this 
sense, the big-bang singularity is already resolved quantum mechanically (see also 241]). 



7. DENSITIZED CONSTRAINT 



Once the state \u = 0) has been removed, let us call Cyl^j the linear span of the nonzero 
volume eigenstates {\u) ; v ^ 0, v G M}. Based on previous experience with gravitational 
models, we expect nontrivial solutions to the constraint to live in the algebraic dual of 
CyVg. Since the operator [l/\/jp[] is invertible in the orthogonal complement of \v = 0), 



it is easy to check that one gets a bijection between the considered solutions and those of 
the alternative "densitized" constraint [3] 



C := -6ft 2 + P 2 . 

The operator ft 2 (with domain Cylg) has the following action: 

fl 2 \u) = - f + (v)f+(v + 2) \v + 4) + [f + {y) + f»] |„) 
- /_(!/)/_(!/- 2) 4), 



(27) 



:28) 



where 



and 



4^ A 1 / 4 



9^) 



g{p±2)s±{v)g{v) 



s±(v) = sign(z/ ± 2) + sign(z/), (29) 



1 

1 + - 

v 



1 

1 - - 

V 



if v ± 0, 



(30) 



while g(z/ = 0) = 0. We notice that VL 2 relates only states \v) whose label differ by a 
multiple of four. Moreover, owing to the linear combination of signs in s±(u), one can see 
that the real function f + (y)f + (y + 2) has a remarkable property, namely, it vanishes in 
the whole interval [-4,0]. Something similar happens with — 2), which vanishes 

in [0,4]. As a consequence, for the label u, the action of the operator ft 2 does not mix any 
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of the semilattices Cf := {±(e + 4n), n G IN}, with £ 6 (0,4] -but otherwise unspecified. 
In the following, we call H.f the corresponding Hilbert subspaces of states with support in 
these semilattices (i.e., the completion of the linear span of ^-states with v e Cf). Each 
of these subspaces can be considered a superselection sector for the quantum theory, 



inasmuch as t 
of the model 



rev provide irreducible representations for the physically relevant operators 
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On the other hand, it is possible to prove that (up to a global multiplicative factor) Q 2 , 
restricted to H+UHj_ e , differs by a symmetric, trace-class operator from an operator which 
is unitarily related with the Hamiltonian of a point particle in a Psch- Teller potential 
22J, [25]. From the properties of this Hamiltonian and Kato perturbation theory 26|], it 
then follows that £7 2 is essentially self-adjoint and that its absolutely continuous spectrum 4 
is 1R + . The rest of the spectrum can be proven empty 22]. Moreover, Vt 2 commutes with 
the projections to H+ and Hj_ £ . One can then see that the operator on any of these 
Hilbert spaces is positive with an absolutely continuous spectrum of unit degeneracy [22 ]. 

Hence, on any superselection sector H^, one obtains a spectral decomposition of the 
identity of the form 

/■oo 

l ±e = / d\\ef*)(e±% (31) 
Jo 

where \e^ e ) is a generalized eigenstate of fi 2 with eigenvalue equal to A. Finally, let us 
comment that, expressing |e^ e ) in the z/-basis, the corresponding generalized eigenfunc- 
tions e^ e (z/) can always be chosen real. 



8. PHYSICAL STATES 

Employing the above spectral decomposition associated with f2 2 , elements of the poly- 
mer space H^ can be identified with elements of the Hilbert space L 2 (M + , dX). It is now 
straightforward to solve the densitized constraint — 6f2 2 + P 2 = 0. Starting from the 
kinematical Hilbert space H^ Cg> L 2 (M, dip), the solutions adopt the form 

poo 

if)(y,<p)= / d\ef{p) V+(A)e iv/ ^ + ?MA)e- iv/ ^ • (32) 



Physical states can be identified with positive frequency solutions, and hence with 
wavefunctions in L 2 (M + ,d\). A complete set of Dirac observables (acting on physical 



4 See, e.g., Appendix C in reference [18] for a summary on operator theory. 
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states) is given by P^ and, e.g., \v\^ , the latter being defined by the action of \v\ when 
= 0o- 5 in this way, the LQC approach succeeds in achieving a complete quantization 
of the flat FRW model with massless scalar field. 

Rather than in general physical states, one is usually interested in states which display 
a semiclassical behavior in the region of large spatial volumes and matter fields, so that 
they can be regarded as potential candidates to explain the properties of universes like 
the one which we observe. With this motivation, we can concentrate our considerations 
on positive frequency states which, for a fixed large value of the scalar field <p — (fio ^> 1, 
are peaked on certain values P$ = P° and v = u° of the Dirac observables such that 
^> 1 and \P®\ 1 [161 ] . In more detail, one analyzes Gaussians of the form 

^ U = ^pj oc e -(-+^°) 2 /( 2 - 2 ) e -^, with 0! = 0o - l/|ln (33) 

Numerical integration of the quantum evolution dictated by the densitized Hamiltonian 
constraint shows that the state remains peaked on a trajectory which coincides with the 
union of a contracting and a expanding classical solutions except in the region where 
the matter energy density becomes comparable to the Planck density in order of magni- 

I 

tude [16j . At that moment, the effective trajectory on which the state is concentrated 
passes from a contraction to an expansion phase in such a way that the classical singu- 
larity is avoided. This quantum phenomenon which allows the resolution of the big bang 
singularity is usually called big bounce [15|, [16J . 



9. CONCLUSION 



We have seen that the quantization techniques of LQG prove to be successful in achiev- 
ing a rigorous and complete quantum theory of simple cosmological models, like e.g. the 
case of flat FRW spacetimes provided with a minimally coupled scalar field. The result- 
ing quantization adopted in LQC is inequivalent to the standard Schrdinger (or Wheeler- 
DeWitt) quantization which has been traditionally employed in geometrodynamics, a fact 
that explains why the emerging physics is radically different and allows LQC to supply 

5 Recalling expression (|32|) . this suffices to determine the action of the operator on positive frequency 
solutions for all values of (p. 
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satisfactory answers to fundamental problems that had remained open in Quantum Cos- 
mology. In particular, this explains why, while the standard quantization fails to solve the 
cosmological singularities, these are cured in LQC. Actually, the singularities are resolved 
already at the kinematical level. Nonetheless, the resolution is much stronger. For physi- 
cal states with good semiclassical behavior, numerical simulations show that the universe 
suffers a big bounce before reaching the big bang. This bounce occurs when the energy 
density p = P|/(2|p| 3 ) approaches a critical density of the order of the Planck density. 
Away from the bounce, states are peaked on classical solutions. Quantum corrections are 
strong close to the bounce, but even there the state remains peaked on a certain, modified 
trajectory. 
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